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We study static properties of disk-shaped binary dipolar Bose-Einstein condensates of 168 Er- 164 Dy 
and 52 Cr- 164 Dy mixtures under the action of inter- and intra-species contact and dipolar interactions 
and demonstrate the effect of dipolar interaction using the mean-field approach. Throughout this 
study we use realistic values of inter- and intra-species dipolar interactions and the intra-species 
scattering lengths and consider the inter-species scattering length as a parameter. The stability of 
the binary mixture is illustrated through phase plots involving number of atoms of the species. The 
binary system always becomes unstable as the number of atoms increases beyond a certain limit. As 
the inter-species scattering length increases corresponding to more repulsion, an overlapping mixed 
state of the two species changes to a separated demixed configuration. During transition from a 
mixed to a demixed configuration as the inter-species scattering length is increased for parameters 
just below the stability line, the binary condensate shows special structures in density in the form of 
red-blood-cell-like biconcave and Saturn-ring-like shapes, which are direct manifestations of dipolar 
interaction. 
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I. INTRODUCTION 

After the experimental realization [THS] of a dipo- 
lar Bose-Einstein condensate (BEC) of 52 Cr atoms with 
magnetic moments, there has been renewed interest in 
the study of static and dynamic properties of such a con- 
densate in the pursuit of novel and interesting properties 
and features emerging as a consequence of anisotropic 
long-range dipolar interaction. The atomic interaction 
in a dilute BEC of alkali-metal atoms is taken as a S- 
wave short-range (delta-function) potential. However, 
the anisotropic long-range dipolar interaction is nonlocal 
in nature acting in all partial waves. More recently, BEC 
of 164 Dy [fH [7] and 168 Er [8] atoms with larger dipole mo- 
ments are available for experimental studies, and polar 
molecules with much larger (electric) dipole moment are 
being considered [3] for BEC experiments. Among the 
novel features of a BEC with anisotropic dipolar inter- 
action, one can mention the peculiar shape and stability 
properties of a stationary state [TU] , red-blood-cell-like bi- 
concave shape due to radial and angular roton-like excita- 
tions [H], anisotropic D-wave collapse [12], formation of 
anisotropic soliton, vortex soliton [T3] and vortex lattice 
[Hj , anisotropic sound and shock wave propagation [TS] , 
and anisotropic Landau critical velocity |16j among oth- 
ers. Distinct stable checkerboard, stripe, and star config- 
urations in dipolar BECs have been identified in a two- 
dimensional (2D) optical lattice as stable Mott insulator 
[17] as well as superfluid soliton [18] states. A new pos- 
sibility of studying universal properties of dipolar BECs 
for large scattering length has been suggested [T!5] ■ Many 
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of these features have only been predicted theoretically, 
although some of them have already been experimentally 
confirmed. To enhance the effect of the anisotropic dipo- 
lar interaction in most of these theoretical studies, the 
contact interaction has been set equal to zero. In 52 Cr 
atom, this could be necessary as the strength of the re- 
pulsive contact interaction is much stronger than that of 
the dipolar interaction. However, dipolar BEC of atoms 
such as Dy with much larger dipolar interaction can 
show effects of anisotropic nonlocal dipolar interaction 
without switching off the contact interaction. 

Now with available dipolar BECs of different species 
of atoms, there is the possibility to investigate a binary 
mixture with two dipolar BECs of two types of dipo- 
lar interactions acting on each species, e.g., intra- and 
inter-species, superposed on intra- and inter-species con- 
tact interactions, thus creating a much richer platform 
to study the effect of dipolar interaction. Previously, 
there have been studies of nondipolar binary boson-boson 
[201121) . boson- fermion [22] and fermion-fermion |23) mix- 
tures. For a binary system without dipolar interactions 
it is possible to have either a mixed or demixed phase. In 
this paper, we consider binary dipolar BECs composed of 
164 Dy and 52 Cr atoms as well as 164 Dy and 168 Er atoms. 
With all the dipole moments aligned along the z axis, 
the net dipolar interaction is attractive in the cigar shape 
along z axis and repulsive in the disk shape confined in 
the x — y plane, as parallel dipoles in a chain along z 
axis attract and those confined in the x — y plane repeal 
each other. Consequently, the binary dipolar BEC with 
dipoles polarized along z axis is more stable in the disk 
shape confined in the x — y plane compared to the cigar 
shape along the z axis vulnerable to collapse and we shall 
consider only such a disk-shaped binary dipolar BEC in 
this study. Even such a disk-shaped binary dipolar BEC 
is found to be unstable due to collapse instability for 
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dipolar interaction above a critical value controlled by 
the number of atoms. This is consistent with a similar 
conclusion that a single-component dipolar BEC for any 
trap anisotropy collapses above a critical value of dipolar 
interaction [2U [25] . 

We study the limits of stability of binary dipolar disk- 
shaped BECs formed of 168 Er and 164 Dy atoms and of 
52 Cr and 164 Dy atoms and present the results in terms of 
stability phase plots which should aid in the experimen- 
tal preparation and study of binary dipolar BECs. Of all 
the atomic interactions, the inter- and intra-species dipo- 
lar interactions controlled by the known dipole moments 
of the two species will be considered known. The con- 
tact interactions governed by the approximately known 
intra-species scattering lengths will also be taken as fixed 
parameters in this study. We then study the stability 
of the binary mixture by varying the number of atoms 
in each species and the yet unknown inter-species scat- 
tering length a\2- For small values of a%2, the dipolar 
binary BEC prefers a mixed configuration and with the 
augmentation of the inter-species repulsion the system 
moves into a demixed configuration. During this process 
of mixing-demixing, distinct biconcave red-blood-cell-like 
and Saturn-ring-like shapes are found in the densities 
of the two components just below the line of stability 
of the binary dipolar BEC for a certain value of inter- 
species scattering length. This configuration corresponds 
to saddle-point structures in the 2D densities of the two 
components in the x — z plane. It was demonstrated 
[3] that the dipolar interaction corresponds to a saddle 
structure in 2D responsible for the saddle structure in 
density in the present study on binary dipolar BEC. The 
biconcave structure is a consequence of roton instabil- 
ity near the stability line due to dipolar interaction and 
was studied in detail in single-component dipolar BECs 
[111 I26j . The Saturn-ring-like density distribution was 
never observed in a single-component dipolar BEC. To 
enhance the dipolar effect, in these previous studies the 
repulsive contact interaction was set to zero, whereas in 
the present study we set all the intra- and inter-species 
repulsive contact interactions to their large experimental 
values and yet obtain these special structures because of 
the larger values of dipole interactions appropriate for 
164 Dy and 168 Er atoms. 

In Sec. II we present the mean-field model for the bi- 
nary dipolar BEC interacting via inter- and intra-species 
contact and dipolar interactions in an axially-symmctric 
confinement. Details of numerics together with an expo- 
sition of our numerical results are included in Sec. III. 
After the application of our model to a simplified binary 
dipolar BEC without any contact interaction, we present 
the results for the binary 168 Er- 164 Dy and 52 Cr- 164 Dy 
mixtures. First, we obtain the stability phase plots of 
the binary dipolar BECs. Then we demonstrate mixing, 
demixing, and structure formation in densities in both 
cases. Close to the stability line, we find Saturn-ring- 
and blood-cell-like density distributions in the compo- 
nents — only possible in the presence of dipolar interac- 



tion. Finally, in Sec. Ill we present a brief summary of 
our findings. 



II. MEAN-FIELD MODEL FOR THE BINARY 
DIPOLAR BEC 

We consider a general two-species dipolar BEC with 
inter- and intra-species dipolar and contact interactions 
with the mass, number, magnetic moment, scattering 
length denoted by m^, Afj, a^, i = 1,2, respectively. 
The angular frequencies for the axially-symmctric trap 
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along x, y and z directions are taken as ui x 
and tu z = AjWj. The intra- and inter-species interactions 
for two atoms at positions ri and r 2 are taken as 
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where R = ri — r 2 , po is the permeability of free space, 9 
is the angle made by the vector R with the polarization 
z direction, and — mimj/fmi + 7712) is the reduced 
mass of the two species of atoms. With these interac- 
tions, the coupled Gross-Pitaevskii (GP) equations for 
the binary dipolar BEC can be written as [1] 
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with p 2 — x 2 + y 2 . 

To compare the dipolar and contact interactions, the 
intra- and inter-species dipolar interactions will be ex- 
pressed in terms the length scales a d l \ and respec- 
tively, defined by 
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For the intra-species dipolar scale aA, the mass of the 
corresponding species rrii has been used to define the 

( 12) 

scale. To define the inter-species dipolar scale a dd , we 
have used the mass toi of the first species. We express 
the strengths of the dipolar interactions in Eqs. ([3| and 
Q by these length scales and transform these equations 
into the following dimensionless form: 
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3JV 2 ai 2) mi2, 5 ^ 2) = 3iV 2 a^ 2) , ^ 21) = SiV^f. In Eqs. 
ufj and ([8]), length is expressed in units of oscillator 
length for the first species Iq = yjh/miusi, energy in 
units of oscillator energy Huii, density |0j| 2 in units of 
Iq 3 , and time in units of to = ^i 1 ■ 



III. NUMERICAL RESULTS 

For the binary dipolar BEC we solve Eqs. |7| and Q 
numerically after discretization [27]. The divergence of 
the dipolar term at short distances has been handled by 
treating this term in momentum (k) space. The dipolar 
integral in (Fourier) momentum space is tackled by the 
following convolution integral [TJ [28] 



J dv'V dd {K)n{v') = J j^e- lk r V dd (k)n(k), (9) 

with n(r) = |0(r)| 2 . The Fourier transformation (FT) is 
defined by 
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The FT V dd (k) of the dipolar potential is known analyt- 
ically [2 [28] : 
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The FT n(k) of density is calculated numerically by a 
fast FT (FFT) routine. The inverse FT is also evaluated 
numerically by the FFT routine. The whole procedure 
is performed in three-dimensional (3D) Cartesian coordi- 
nate system irrespective of the underlying trap symme- 
try. 



A. Same species of atoms in isotropic trap 

To test the model and the numerical routine, we 
consider a simple binary dipolar BEC of same species 
of atoms without contact interaction in a spherically- 
symmetric trap. However, the two components are con- 
sidered to be oppositely polarized along z and — z di- 
rections, respectively. In this model, first considered by 
Goral and Santos 

a dd ,ai = 012 = 0, tot 2 = niuj = 1,A, ; = 1. Con 
scqucntly, in Eqs. 
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3iVjOdd, g dd "' = -3N 2 a dd , and g dd = -3Nia dd . Using 
r] = Ni/N, C = 3a dd N, N = N 1 + N 2 we can write Eqs. 
and (|8| as in Ref. [25]: 
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For a fixed fraction of the atoms of the first component 
77, the system is stable up to a critical value of net dipo- 
lar interaction nonlinearity Cj e.g. for £ < £ cr , beyond 
which the system becomes unstable. The stability of the 
binary system can be expressed in terms of a C — V phase 
diagram. This stability phase diagram for this binary 
dipolar condensate of same species of atoms, but of op- 
positely polarized dipolar moments, in a spherically sym- 
metric trap is shown in Fig. 1 in agreement with Goral 
and Santos [25J . Independent of the composition of the 
system governed by the parameter 77, the system is unsta- 
ble beyond a critical value of the total number of atoms 
N. 



B. Different species of atoms in disk trap 

Binary 168 Er- 164 Dy mixture: We consider two dif- 
ferent binary mixtures of dipolar BECs. First, we con- 
sider the 168 Er- 164 Dy mixture. In this case we take BEC 
number 1 as 168 Er and BEC number 2 as 164 Dy with 

parameters [i\ = 7/is,/i 2 = 10/is, a^J = 660,0,0,^^ = 
(2) 

94ao,a dd = 131ao with ao the Bohr radius and [Ib the 
Bohr magneton. Without accurate experimental esti- 
mates [7[ [5J of the intra-species scattering lengths, we 
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FIG. 1: (Color online) ( — r\ stability plot for a binary con- 
densate of oppositely polarized dipolar gases in a spherically 
symmetric trap compared with the findings of Goral and San- 
tos EE- 



use di = llOao. The angular frequencies of the axial 
trap for 164 Dy are taken as u>2 — 2-7T x 243 Hz, A2 = 
\/l0 ~ 3.1623, compared with experimental frequencies 
H {fx, fy, fz} = {205, 195, 760} Hz with trap aspect ra- 
tio A = 3.8. We take m w = (m 2 w|)/(miw 2 ) = 1, corre- 
sponding to lji — 2-7T x 240 Hz, Ai = %/TO. The unit of 
length in this study is Iq — \Jhj(m\ijj\) rj 0.5 /xm and 
the unit of time t = uj^ 1 w 0.663 ms. 

Binary 52 Cr- 164 Dy mixture: In this case, taking the 
condensate number 1 as 52 Cr and the one number 2 
as Dy, the parameters are: [i\ = 6 \i b, H2 — 10/i_b, 



l dd 



= 15a o,ald 
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= 25ao- We take ai = 



110a , close to their experimental estimates H |7j. 
The angular frequencies for the axial trap for 52 Cr are 
taken as w\ = 2tt x 195 Hz, Ai = vlO « 3.1623, com- 
pared with the experimental frequencies [4] {fx, fy, fz} — 
{150,150^/^} Hz with trap aspect ratio A < 10. We 
take m u = (m2u; 2 )/(miu; 2 ) = 1, corresponding to cj 2 — 
2tt x 110 Hz, A2 = \/T0. The unit of length in this calcu- 
lation is Iq = \J h/{m\(jj{) ~ 1 /im and the unit of time 
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C. Stability phase plot 

By solving the coupled set of GP equations ([7 ) and Q 
for the binary mixture of dipolar BECs with the above 
parameters we find that the system becomes unstable 
above a certain number of atoms with the total dipolar 
interaction beyond a limiting value in agreement with 
a similar conclusion [24] in a single-component dipolar 
BEC. To perform a systematic study of the stability of 
the binary dipolar BEC, we solve Eqs. (JtJ) and (JsJ) with 
the above parameters, fixing the number of atoms in the 
first species ( 52 Cr or 168 Er) and searching for the critical 
number of atoms (N cr ) of second species 164 Dy beyond 
which the system becomes unstable for different values of 
inter-species scattering length ai2. If the binary mixture 
has a smaller number of atoms of the first species ( 52 Cr 
or 168 Er), it can accommodate a larger number of atoms 
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FIG. 2: (Color online) Stability phase plot showing the critical 
number N cr (Dy) of 164 Dy atoms in the 52 Cr- 164 Dy and 168 Er- 
164 Dy binary mixtures versus inter-species scattering length 
ai2 for (a) 10000 and (b) 1000 52 Cr or 168 Er atoms, respec- 
tively. The system is stable below the respective lines. The 
shaded dark areas illustrate the metastable domains, where 
biconcave density profiles appear. 



of the second species ( 164 Dy) while the limiting value of 
the total dipolar interaction is reached maintaining all 
other parameters fixed. 

We present the results of our study in stability phase 
plots in Fig. 2. In this figure we show the critical num- 
ber of 164 Dy atoms N cr (Dy) versus a u for 100 52 Cr or 
168 Er atoms in binary 52 Cr- 164 Dy and 168 Er- 164 Dy mix- 
tures, respectively. The system is stable below the lines 
of Fig. 2 and unstable above. Stability phase plots of 
Fig. 2 could be relevant for planning future experiments 
on binary dipolar BEC. The system is only metastable in 
the shaded regions in this figure, possibly corresponding 
to a local and not a global minimum of energy. This re- 
gion is similar to the darker region in Fig. 1 of Ref. |24j . 
where biconcave shape in density of a single-component 
dipolar BEC appeared due to roton instability. The 
biconcave shape in density in a disk-shaped dipolar BEC 
emerges because of dipolar repulsion in the plane of the 
disk. Due to this repulsion the atoms come to the periph- 
eral region of the disk and a region of low density appears 
in the center. Similar low-density central region may ap- 
pear in a rotating BEC due to centrifugal repulsion, as 
has been found in a binary nondipolar BEC in a rotat- 
ing trap [H]. We shall see below that biconcave and 
Saturn-ring-like shapes in density of the components of 
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the binary dipolar BEC may appear in the shaded regions 
in Fig. 2. The system can accommodate a larger number 
of atoms for large positive values of an responsible for 
large contact repulsion which stabilize the binary dipo- 
lar BEC. As a%2 reduces, the contact repulsion reduces 
and the system becomes more vulnerable to collapse for 
larger number of atoms due to dipolar interaction and 
hence can accommodate only a small number of 164 Dy 
atoms as can be seen in Fig. 2. 



D. Mixing, demixing, and structure formation 

Next we study in detail the density of the two species 
of atoms in the binary dipolar BEC to look for mixing 
(overlapping phases of two components) , demixing (sepa- 
rated phases) and structure formation in the metastable 
regions of the phase plots shown in Fig. 2. For small 
ai2 the inter-atomic contact repulsion between the two 
species is small and one has a mixed configuration. As a\2 
is increased the inter-species contact repulsion increases 
and for a large enough value of a\% one can have demix- 
ing of the two species. In the demixed configuration one 
of the species occupy the central region and the other 
the external region. To illustrate this phenomenon, we 
plot in Fig. 3 the reduced 2D density in the radial direc- 
tion \$(x,y = 0)| 2 — J dz\4>{x,y — 0, z)\ 2 versus x under 
different situations in case of the binary 168 Er- 164 Dy mix- 
ture. The mixing-demixing is illustrated in Figs. 3 (a), 

(b) , and (c) for 10000 168 Er and 20 00 164 Dy atoms, re- 
spectively, for ai2 = llOao, 125ao, and 130ao. In Fig. 3 

(a) we have a mixed (overlapping) configuration and with 
the increase of inter-species contact repulsion, in Fig. 3 

(c) we have a demixed (separated) configuration of the 
two species. In Fig. 3(b), for a 12 — 125ao we have a 
partially demixed configuration. In the demixed configu- 
ration in Fig. 3 (c), 20 00 164 Dy atoms have a larger con- 
tribution to energy of the system and the 164 Dy atoms 
expel the 168 Er atoms from the central region due to 
inter-species repulsion. The opposite panorama is also 
possible for a larger number of 168 Er atoms, while the 
168 Er atoms can expel the 164 Dy atoms from the central 
region. Such a situation is illustrated in Fig. 3 (d) with 
100 00 168 Er atoms and 500 164 Dy atoms for a 12 = 160a Q . 
In this case 168 Er atoms have a larger contribution to the 
energy of the system and expel the 164 Dy atoms from the 
center. 

A careful examination of the passage from the partially 
mixed configuration of Fig. 3 (b) to the demixed config- 
uration of Fig. 3 (c) reveal interesting feature. To study 
this, we consider the 3D profile of the density \4>i(x, y, z)\ 2 
of the two species for the cases depicted in Figs. 3 (b) 
and (c). In Figs. 4 (a) and (b) we plot the 3D contours 
of 168 Er and 164 Dy BEC profiles for the parameters of 
Fig. 3 (b). The same for the parameters of Fig. 3 (c) 
are illustrated in Figs. 4 (c) and (d). In Figs. 4 (a) and 

(b) the 168 Er and 164 Dy BECs have Saturn-ring-like and 
red-blood-cell-like biconcave profiles, respectively. The 
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FIG. 3: (Color online) Two-dimensional radial density along 
x axis \&(x,y = 0)| 2 = J dz\(fi(x,y — 0,jz)| 2 of the binary 
i68 Er _i64 Dy BEC] for ^ N ^ = 1000Q] 7v(Dy) = 20000, 

and Oia = 110a ; (b) N(Er) = 10000, N(Dy) = 20000, and 
ai2 = 125a ; (c) AT(Er) = 1000, JV(Dy) = 20000, and a 12 = 
130a ; (d) iV(Er) = 10000, JV(Dy) = 500, and a 12 = 160a . 
For the trap parameters of this study the length scale lo = 0.5 
fj,m. 



biconcave profile is a manifestation of the dipolar inter- 
action: the net dipolar repulsion in 164 Dy in disk shape 
drives the dysprosium atoms to peripheral region thus 
creating the biconcave shape as in the single-component 
dipolar BEC. In the single-component case the bicon- 
cave profile near the stability limit has been related to 
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FIG. 4: (Color online) Three-dimensional contour plot of the 
(a) 168 Er and (b) 164 Dy densities \(f>i(x, y, z)\ 2 of the 168 Er- 
164 Dy binary BEC of Fig. 3 (b) for a 12 = 125a . The same 
for the (c) 168 Er and (d) 164 Dy densities of the binary BEC 
of Fig. 3 (c) for a\2 = 130<zo. Density on contour = 0.0005, 
and 7V(Er) = 10000, JV(Dy) = 20000. The lengths x, y, z are 
in units of Iq = 0.5 /J,m. 



roton instability [24]. The Saturn-ring- like 168 Er profile 
in Fig. 4 (a) is a manifestation of the inter-species dipo- 
lar interaction: the biconcave density distribution of the 
164 Dy species with a density minimum at center, due to 
inter-species repulsion, expels a part of the 168 Er BEC to 
the central region of lower density and another part to 
the peripheral region, also of lower density, thus creating 
the Saturn-ring-like profile. Such a profile is not possi- 
ble in a single-component dipolar BEC, or in a binary 
BEC without dipolar interaction. With further increase 
of inter-species contact repulsion, for an — 130ao, the 
contact repulsions dominate over the dipolar interactions 
and hence plays a major role thus creating a simpler sit- 
uation of demixing shown in Fig. 3 (c) and Figs. 4 (c) 
and (d) , where the net inter-species repulsion is so strong 
that all 168 Er atoms are driven to the peripheral region 
in the form of a ring with the dominant 164 Dy atoms 
occupying the central region. The ring-like structure in 
Fig. 4 (a) as opposed to a shell-like demixed configura- 
tion is a consequence of the dipolar interaction: a pure 
inter-species contact repulsion would have led to a hollow 
shell-like configuration for 168 Er surrounding a compact 
disk-shaped 164 Dy core. The dipolar interaction is attrac- 



tive in the z direction and this transforms the shell-like 
configuration into a ring. Such transition from a shell- 
like configuration to a ring with the increase of dipolar 
interaction has been demonstrated in a single-component 
dipolar BEC in a shell-like trap [29] . 

From Figs. 3 we find that for a binary dipolar BEC 
composed of 10000 168 Er atoms and 20 00 164 Dy atoms 
the transition from a mixed to a demixed configuration 
happens around a 12 = 120a . We verified that simi- 
lar transition also takes place in a binary dipolar BEC 
for 100 00 168 Er atoms and a reduced number of 164 Dy 
atoms also for a\2 ~ 120ao- Nevertheless, the biconcave 
structure in the density of 164 Dy is obtained for a small 
shaded region below the stability line in Fig. 2 (a) with 
large number of 164 Dy atoms. As in the single component 
case [24], the the binary dipolar BEC possessing bicon- 
cave shape in density in the 164 Dy BEC in the shaded 
region in Fig. 2 (a) is at best metastable corresponding 
to a local minimum (and possibly not a global minimum) 
in energy and instability appears only after crossing the 
stability lines in Fig. 2. As pointed out in the case of 
a single-component dipolar BEC a conveniently chosen 
initial density could be necessary for obtaining the bicon- 
cave structure in density. In imaginary-time propagation 
the initial density should either be chosen as a Gaussian 
with widths slightly larger than the actual widths of the 
desired state [30] or as random initial functions [26]. In 
case of random initial function one needed to take an an 
average over 100 realizations of imaginary time propaga- 
tion to obtain the biconcave structure. The present struc- 
tures were obtained by using imaginary time propagation 
of the coupled mean-field equations with Gaussian pro- 
files for initial densities. We started from the Gaussian 
profiles of the linear oscillator states and increased the 
dipolar and nondipolar nonlinearities in small steps dur- 
ing time evolution until the final nonlinearity is reached. 
In this fashion the metastable states in the shaded regions 
of phase plots in Fig. 2 are obtained. As these structures 
are only metastable, they may disappear upon evolution 
for a very long time. 

To gain further insight into the density distributions in 
Figs. 4 (a) — (d) we plot in Figs. 5 (a) — (d) contour plots 
of asymmetric 2D density |$(a;,z)| 2 = J dy\(f>(x, y, z)\ 
of the binary 168 Er- 164 Dy mixture for parameters corre- 
sponding to Figs. 4 (a) — (d), respectively. The double 
and single saddle structures in density in the x — z plane 
in Figs. 5 (a) and (b) are manifestations of the dipolar 
interaction. In these plots a minimum in density along 
x direction coincides with a maximum in density in z di- 
rection thus creating a saddle point, which appears as a 
clear manifestation of the saddle-shaped dipolar interac- 
tion [3]. Such a density distribution is not possible in 
the absence of dipolar interaction and to the best of our 
knowledge has not been demonstrated even in the single- 
component dipolar BEC. In Fig. 5 (a), for a\% = 125ao, 
we have two saddle points on the x axis in the density of 
168 Er and one saddle point in the density of 164 Dy. Simi- 
larly, for a\2 — 130ao, we have a single saddle point only 
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FIG. 5: (Color online) Contour plot of 2D density z)| 2 = 
Jdi/|^(x,y,2)i 2 of the binary 168 Er- 164 Dy BEC for JV(Er) = 
10000, iV(Dy) = 20000 with a 12 = 125a for (a) 168 Er, (b) 
164 Dy, and with a 12 = 130a for (c) 168 Er, and (d) 164 Dy. 
The lengths are expressed in units of Iq(= 0.5 fim) and the 
2D density in units of Iq 2 . 



in the density of 168 Er and none in 164 Dy. These sad- 
dle points have curvatures consistent with saddle-shaped 
dipolar interaction [3] with negative curvature along the 
direction of magnetization (z axis) and positive curvature 
orthogonal to it (x axis). 

In Fig. 6 we plot the chemical potential fi and root 
mean square (rms) sizes (z), and (p) versus inter-species 
scattering length o 12 , for (a) JV(Er) = iV(Dy) = 10000; 
(b) iV(Er) = 10000, iV(Dy) =1000; (c) iV(Er) = 1000, 
N(Dy) = 10000. It is found that the chemical potential 
and rms sizes are slowly varying functions of inter-species 
scattering length. As expected, the chemical potential 
of the first species increases more rapidly with a\2 when 
there are more atoms in the second species and vice versa. 

With this study of the binary 168 Er- 164 Dy BEC, we 
also considered a detailed investigation of the binary 
52 Cr- 164 Dy BEC. Similar mixing, demixing, and struc- 
ture formation are also found in this second case, al- 
though the quantitative results and estimates are differ- 
ent. We only report here the studies on the interesting 
structure formation in the binary 52 Cr- 164 Dy BEC. In- 
deed we find that one can have the condensate with bi- 
concave shape in this case. In Fig. 7 (a) and (b) we plot 
the linear density in the x — y plane along the x direc- 
tion \$>{x,y = 0)| 2 , as in Figs. 3, for a binary mixture 
with 100 52 Cr atoms and 60 00 164 Dy atoms and for 
ax2 = 70a and 80ao, respectively. Figure 7 (a) illustrates 
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a case similar to Fig. 3 (b) denoting a transition from 
a mixed state to a demixed state. For a\2 < 60ao, we 
have a fully mixed configuration of the two condensates 
as in Fig. 3 (a). For = 70ao a biconcave shape has 
appeared in the density of 52 Cr. For a\i > 75ao, a com- 
plete demixed configuration appears as illustrated in Fig. 
7 (b) for a = 80ao, where the 52 Cr atoms are expelled 
from the central region occupied by only 16 Dy atoms. 
In Figs. 7 (c) and (d) we show the 3D contour of the 
52 Cr and 164 Dy BECs corresponding to the parameters 
of Fig. 7 (a). The 52 Cr BEC shows a pronounced bicon- 
cave shape which is a direct manifestation of the dipo- 
lar interaction. However, we could not find a biconcave 
shape in the density of the 164 Dy BEC, as in the case of 
a binary dipolar 168 Er- 164 Dy mixture, possibly because 
the inter-species dipole interaction in the present system 
is much smaller in strength compared to the same in the 
binary 168 Er- 164 Dy system so as to create the biconcave 
shape in density in the presence of dominating large con- 
tact repulsions acting on 164 Dy. The much larger dipole 
interactions in 164 Dy might be responsible for the bicon- 
cave shape in density in Fig. 4 (b). Even in the single- 
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FIG. 7: (Color online) Two-dimensional radial density along 
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= 50000, for (a) 



and (c) Dy density \4>%(x,y 
Density on 3D contour = 0.0025 
units of l = 0.5 fim. 



2 of the binary BEC of (a). 
The lengths x, y, z are in 



of 164 Dy BEC is shown in Fig. 8 (c). The central high- 
density region of 164 Dy BEC has expelled the 168 Er BEC 
from the central region, thus creating a biconcave shape 
in density. If we compare Figs. 4 (a) and (b) with Figs. 
8 (b) and (c), we find that the roles of 164 Dy and 168 Er 
BECs have interchanged. In Fig. 4 (a), the 168 Er BEC 
has a Saturn-anel-like profile, whereas, in Fig. 8 (c), 
the 164 Dy BEC has similar profile. These profiles are a 
consequence of dipolar interaction. 



component case, the biconcave profile in density appears 
for certain values of the trap aspect ratio [24], and it is 
also possible that for other values of trap aspect ratio the 
biconcave profile might appear in the density of 164 Dy in 
the present binary 52 Cr- 164 Dy mixture. 

Finally, we study the metastable binary dipolar BECs 
with a smaller number 168 Er atoms in the shaded region 
of Fig. 2 (b). In particular we consider the partially 
mixed configuration with an — 110ao,7V(Er) = 1000, 
7V(Dy) = 50000. In Fig. 8 (a), we illustrate the 2D radial 
density along x axis y — 0)| 2 = J dz\<f)(x, y = 0, z)\ 2 
of this binary 168 Er- 164 Dy BEC. In Figs. 8 (b) and (c) 
the contour plots of the corresponding 3D densities of 
168 Er and 16 Dy are shown. In this case a biconcave 
shape has appeared in the density profile of 168 Er BEC 
in Fig. 8 (b). A Saturn-anel-like profile in the density 



IV. SUMMARY 

Using a mean-field description we studied the static 
properties of binary disk-shaped dipolar BECs 168 Er- 
164 Dy and 52 Cr- 164 Dy to search for the effect of inter- 
and intra-species dipolar interactions employing realis- 
tic values of inter- and intra-species dipolar interactions 
as of intra-species scattering lengths. The yet unknown 
inter-species scattering length is considered as a variable 
parameter. The binary system is found to be stable for 
number of atoms below a critical value. The stability 
domain for the system is illustrated in convenient phase 
plots involving number of atoms of the species and the 
inter-species scattering length a\i. For small values of 
an a mixed state (with overlapping phase) of the two 
species emerges and this state transforms into a demixed 
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state (with separated phase) for larger a 12 . Just be- 
low the stability line, in the region of transition from 
the mixed to a demixed configuration special metastable 
structures in 3D densities may appear in the form of 
the Saturn-ring-like shape or red-blood-cell-like bicon- 
cave shape corresponding to a local minimum in energy. 
Such structures are a direct manifestations of the dipo- 
lar interaction. Similar biconcave shape in the density of 
a single-component dipolar BEC was found and related 
to roton instability, which is presumably also responsible 



for these structures in the binary dipolar BECs. After 
this investigation was finished we came to know about a 
complimentary work on binary dipolar BEC in 2D |31j . 
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